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ABSTRACT 

We present the first results of tlie VLBA Imaging and Polarimetry Survey 
(VIPS), a 5 GHz VLBI survey of 1,127 sources with flat radio spectra. Through 
automated data reduction and imaging routines, we have produced publicly avail- 
able I, Q, and U images and have detected polarized flux density from 37% of 
the sources. We have also developed an algorithm to use each source's I image to 
automatically classify it as a point-like source, a core-jet, a compact symmetric 
object (CSO) candidate, or a complex source. Using data from the Sloan Digital 
Sky Survey (SDSS), we have found no significant trend between optical fiux and 
5 GHz fiux density for any of the source categories. Using the velocity width of 
the H/5 emission line and the monochromatic luminosity at 5100 A to estimate 
the central black hole mass, Mbh, we have found a weak trend between Mbh and 
5 GHz luminosity density for objects with SDSS spectra. Ongoing optical follow- 
up for all VIPS sources will allow for more detailed explorations of these issues. 
The mean ratio of the polarized to total 5 GHz fiux density for VIPS sources with 
detected polarized fiux density ranges from 1% to 20% with a median value of 
about 5%. This ratio is a factor of ~3 larger if only the jet components of core- 
jet systems are considered and is noticeably higher for relatively large core-jet 
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systems than for other source types, regardless of which components (i.e., core, 
jet, or both) are considered. We have also found significant evidence that the 
directions of the jets in core-jet systems tend to be perpendicular to the electric 
vector position angles (EVPAs). The data is consistent with a scenario in which 
~24% of the polarized core-jets have EVPAs that are anti-aligned with the di- 
rections of their jet components and which have a substantial amount of Faraday 
rotation. Follow-up observations at multiple frequencies will address this issue 
in more detail. In addition to these initial results, plans for future follow-up 
observations are discussed. 

Subject headings: galaxies: active - surveys - catalogs - galaxies: jets - galaxies: 
nuclei - radio continuum: galaxies - techniques: image processing 



1. Introduction 

Very Long Baseline Interferometry (VLBI) can be a powerful tool for the detailed study 
of the nature of the centers of active galaxies. Because VLBI can provide parsec-scale images 
at large distances and because active galactic nuclei (AGN) are among the few objects that 
have brightness temperatures that are high enough to be detected with VLBI instruments. 



VLBI surveys of AGN have been conducted (e.g.. 


Pearson &; Readhead 


1988 


1996; 


Kellermann et al. 


1998: 


Beaslev et al. 


2002: 


Lister & Homan 


2005). 



Taylor et al. 



have pioneered the effort to more fully understand the nature of AGN on parsec scales, but 
still do not provide both imaging and polarization data of samples large enough to answer 
many key questions such as: 



How do the nature and properties of jets evolve as they propagate from their black hole 
sources through different size scales out to the spatial extent of radio lobes (i.e., from 
~10 AU up to ~10 kpc)? Observations of AGN in the X-ray regime with the Chandra 
satellite will help explore this issue on relatively large (~1 kpc) scales. Observations of 
emission within the gamma-ray regime conducted with the High Energy Stereoscopic 
System (HESS) have provided insights into the nature of the gamma-ray emittin g 
regions (~0.01-1 pc) of AGN, but mostly for BL Lac objects (lAharonian et al. II2005I ) . 
A VLBI imaging survey combined with monitoring by the upco ming (launch in 2007) 



Gamma-Ray Large Area Space Telescope (GLAST) mission (iGehrels fc Michelson 



19991 ). and follow-up observations across the electromagnetic spectrum will be able 
to address this issue on relatively small (~10 pc) scales for a variety of AGN. 
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How are synchrotron radiation-emitting particles accelerated along jets and are these 
jets confined by toroidal magnetic fields or gas pressure? A combination of radio spec- 
tral studies, multi-wavelength observations, and high spatial resolution polarimetric 
imaging would help answer these questions. A VLBI survey with imaging and po- 
larimetry for a large sample of AGN is the first key step in this process. 

Is there a statistically significant trend between the direction of core magnetic fields 
and the direction of jets a mong different classes of AGN? While s uch trends have 
been observed for quasars (jPoUack et al. 1 120031 ) and BL Lacs (e.g., iGabuzda et al. 
2OOOI ). the lack of good polarization information for both classes of sources weakens 
the significance of these observed trends. A relatively large sample of AGN with such 
polarization data would help this effort immensely. 



How do radio sources associated with the central black holes of galaxies evolve and af- 
fect galaxy evolution? A key to answering this question may be the subclass of objects 
known as Compact Symmetric Objects (CSOs) which may evolve into sources that 
resemble more well k nown classes of radio galaxies. However, the relatively small num- 
ber of known CSOs (jPeck fc Taylor II2OOOI ) precludes any definite conclusion regarding 
their evolution. A large sample of candidate CSOs imaged using VLBI at multiple 
epochs would help answer this question. 



In order to address these key questions, we have compiled an imaging and polarimetry survey 
of 1,127 AGN with the NRAO Very Long Baseline Array (VLBA). The VLBA Imaging and 
Polarimetry Survey (VIPS) consists of images and polarization data at 5 GHz with follow-up 
observations planned at different epochs and frequencies for selected sources. Several VIPS 
sources will also be found to flare by GLAST, and will subsequently be the target of further 
VLBI follow-up campaigns to connect the presumed jet ejection with the gamma-ray flare. 
This paper describes the sample selection and VLBA observations (§2.1), the data reduction 
and automated imaging (§2.2), and the classiflcation of sources (§2.3). A discussion of flrst 
results regarding the fraction of polarized flux density among the sources and a summary of 
future follow-up plans are contained in §3. 



2. Sample Definition and Data Processing 

2.1. Sample Selection and Observations 

To meet the primary goals of this project, a relatively large sample of likely AGN, 
preferably with data from other wavelength regimes, is required. To this end, we have have 
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chosen the Cosmic Lens All-Sky Survey (CLASS; iMyers et al. II2003I ) as our parent sample. 
CLASS is a VLA survey of ~12,100 flat-spectrum objects (a > —0.5 between 4.85 GHz and 
a lower frequency), making it an ideal source of likely AGN targets to be followed up with 
the VLBA. We have also restricted o ur sample to lie on the survey area, or "footprint" of 
the Sloan Digital Sky Survey (SPSS: lYork et al. 1120001 ). Through the fifth data release of 
the SDSS (DR5; lAdelman-McCarthy et al. II2006I ). the imaging covers 8,000 square degrees 
and includes ~ 2 x 10® objects. Spectroscopy was obtained as part of the SDSS for ~ 10^ 
of these objects, about 10^ of which are quasars. We have chosen our source catalog so that 
all sources lie on the original SDSS footprint with an upper declination limit of 65° imposed 
to avoid the regions not imaged through DR5 (see Fig. [1]). We also excluded sources below 
a declination of 15° because it is difficult to obtain good {u, v) coverage with the VLBA for 
these objects. To keep the sample size large but manageable and to obtain a high detection 
rate without phase referencing, we selected all CLASS sources within this area on the sky 
with flux densities at 8.5 GHz greater than 85 mJy, yielding a sample of 1,127 sources. 
Among these sources, 1,043 (93%) have SDSS images and 356 (32%) have SDSS spectra 
through DR5 with a median redshift of 1.2. 

Among our target list, 141 have already been observe d at 5 GHz with the VLBA as part 



of th e Caltech-Jordell Ban k Flat spectrum survey (CJF; iTaylor et al. Ill996l : iBritzen et al. 



20031 : iPoUack et al. II2003I ). 8 have been or will be observed at 15 GHz as part of the Mon i- 



toring of Jets in AGN with VLBA Experiments project (MOJAVE; iLister fc Homan 
and 20 were observed for the VIPS pilot program at 5 and 15 GHz (jTaylor et al~ 



2005) 



2005|). 



Each of these surveys contain VLBA observations at 5 or 15 GHz that have adequate sensi- 
tivity for our pur poses. Both the MOJA VE and the VIPS pilot surveys were observed in full 
polarization, and iPoUack et al. I (120031 ) obtained full polarization data for over half of the 
sources from the CJF survey that are also within our sample. Because of this, we have elected 
to not re-observe sources from our sample contained within the CJF, MOJAVE, and VIPS 
pilot surveys. The remaining 958 sources were observed with the VLBA within 18 separate 
observing runs of approximately 11 hours each from January to August 2006. The targets 
for the runs consisted of groups of 52-54 VIPS catalog sources with four separate calibration 
sources, 3C279, J1310+3220, and some combination of DA193, OQ208, 3C273, J0854+2006, 
and J1159+2914. Each VIPS target was observed for approximately 500 seconds divided into 
10 separate scans. All observations were conducted with four 8 MHz wide, full polarization 
IFs centered at frequencies of 4609, 4679, 4994, and 5095 MHz. For these observations, an 
aggregate bit rate of 256 Mbps was used, yielding increased baseline sensitivity relative to 
the pilot program observations which used at a bandwidth of 128 Mbps. 

All VLBA observations were scheduled using version 6.05 of the VLBA SCHED program. 
Using built-in data regarding the locations and operation of the VLBA stations, a new mode 
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in the updated version of SCHED will automatically produce a schedule for a list of targets 
with scan durations, a starting LST, and total experiment duration that is optimized both 
for [u, v) coverage and efficiency. For each observing run, the starting LST and scan time 
per source was varied to produce a schedule that most efficiently used the entire duration of 
11 hours while obtaining the vast majority (if not all) of the required scans for all targets. 
Care was also taken to select the correct polarized calibration source(s) for each run so that 
it / they would be observed over a wide range of parallactic angle values while not significantly 
reducing the efficiency of the schedule for that run. 



2.2. Calibration and Source Mapping 



The bulk of the data calibration and imaging was performed in an automated way using 
both AIPS and DIFMAP scripts, similar to the process used in the VIPS pilot program. The 
i nitial calibration and flag ging of bad data was done using the VLBA data calibration pipeline 
( ISjouwerman et al. II2005I ) in AIPS. All of the Stok es I, Q, and U li nages were generated using 



DIFMAP scripts which are described in detail in iTaylor et al. I (120051 ). Polarized intensity 



and polarization angle images were made in AIPS. For sources with data from only the CJF 
survey, the data was obtained from the NRAO archive and images were produced using the 
same procedure. The basic procedure used is as follows: 



1. The initial calibration is done with the VLBA pipeline. 

2. Using this calibration, maps of the four calibration sources are made with the DIFMAP 
scripts. 

3. In AIPS, the calibration is reflned using the maps produced in step 2 as models for the 
calibration sources for self-calibration of the phases, and then again for self-calibration 
of both the phases and the amplitudes of the visibilities. Polarization corrections are 
also determined using either DA193 or OQ208 and are applied to the data. 

4. Using the new calibration, maps are made of all the sources with the DIFMAP scripts. 
These maps are then used to identify the 16 brightest sources (excluding the calibration 
sources), which typically included all sources with peak flux densities ~100 mJy or 
higher. We chose to use the same number of calibration sources for each observing run 
rather than use all sources above a particular flux density limit so that the quality of 
the calibration would remain roughly constant among the runs. Following the addition 
of new cahbration sources, step 3 is then rerun using the maps of these sources for self- 
calibration (effectively yielding 20 calibration sources) to further reflne the calibration 
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solution. This typically eliminates phase errors due to differences in right and left 
polarization and improves the amplitude calibration. 

5. The position angle of 3C279 C4 is used to determine the corrections needed to align the 
observations at the four different frequencies. After applying these corrections, maps 
of the calibration sources are remade with the DIFMAP scripts and the integrated 
polarization angle is measured for each of them using the Q and U maps for each 
of the four frequencies. Observations of these sources from the 5 GHz VLA/VLBA 
polarization calibration database (see the acknowledgments) are then used to determine 
the mean phase correction needed to align the observed polarization angles with the 
true electric vector position angle (EVPA). 

6. After applying the EVPA correction, final images are made of all sources with the 
DIFMAP scripts with a fixed restoring beam for all images. For these maps, the 
minimum dynamic range required by the scripts for peaks to be identified within each 
image is lowered from 6 in previous applications to 5.5. Q and U maps are made in all 
four frequencies as well as in pairs of frequencies, namely 4609 plus 4679 MHz and 4994 
plus 5095 MHz. This was done to improve the signal-to-noise ratio of the polarized flux 
density images (see the following step) by combining observations that are relatively 
close (within 100 MHz) of each other. While not used here for polarization related 
measurements, the Q and U images made using each of the four different frequencies 
have been made publicly available (see below) for future applications. It should be 
noted that the frequencies used do not provide enough separation in for us to be 
able to obtain precise rotation measures (RMs) with which corrections for Faraday 
rotation may be applied. 

7. The final I, Q, and U maps are used within AIPS to make images of both polarized 
flux density and polarization angle for the two pairs of frequencies using the COMP 
task including error biasing. For any object that has polarized flux density that is 
significantly (> 5cr) larger than the noise in its polarized flux density map, a contour 
map of the fraction of polarized flux density and polarization angle is made using the 
4994 plus 5095 MHz pair. It should be noted that no Stokes V images were created for 
our sources since most quasars do not have detectable amounts of circular polarization, 
and in instances where it is detected, the level o f circular polarization is t ypically much 



lower than the level of linear polarization (e.g., iHoman fc War die Ill998l ). 



The typical noise measured from the 5 GHz I images, rmsjmage, is about 0.2 mJy beam ^. 
We have also computed the theoretical noise, rmstheory, according to 

SEFD 

rmStheory = = (1) 
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This equation holds for the simphfied case of Stokes I images generated with N^s visibihties 
measured with identical antennas with natural weighting and no tapering, where SEFD 
is the system temperature in units of Jy, 1]^ is the system efficiency, DR is the data rate 
in bits s~^, and is the time interval in seconds over which the visibi hties were average 



to pr oduce the N^s visibilities that were used to generate the image (IWrobel fc Walker 



9991). For the VL BA, the average value of SEFD for the 10 antennas is about 300 Jy 



(ITaylor et al. Ill994i ). As mentioned earlier, for the newly observed VIPS sources, the data 



rate is 256 Mbps as compared to the data rate of 128 Mbps used for the pilot survey. For 
those sources with data only from the CJF, the data rate ranges from 16 to 64 Mbps and was 
computed for each individual source assuming 1-bit sampling and using the total bandwidth 
and the number of independent polarizations used. For the CJF data, we assumed a system 
efficiency of l/rjs = 1.8 for the case of 1-bit sampling. For the newly observed VIPS sources 
and those sources from th e pilot survey, we assumed l/i]s = 1.84 for the case of 2-bit sampling 
JWrobel fc Walker I Il999h . 

In Fig. [21 we have plotted the ratio of rmsimage to rmsthcory versus the peak signal- 
to-noise ratios for the images. More than 80% of the newly observed VIPS sources have 
rmSimage >i'niStheory and nearly all of these sources have rmSimage < 1-7 rmstheory This 
implies that in general, the newly acquired data is not dynamic range limited and that our 
automated use of the CLEAN algorithm has not artificially reduced the noise level of the 
images. For the data taken from the CJF survey, however, more than half of the sources have 
rmSimage <i'mstheory While our simplified computation of rmstheory has not taken into account 
factors such as weather and elevation, a proper treatment of such effects would most likely 
increase the values of rmstheory and worsen the discrepancy between rmsimage and rmstheory 
We believe that the quality of the CJF data has primarily contributed to the relatively 
low values of rmSimage as the images produced from these data can have considerable peaks 
within the noise which may lead to the images being "over-cleaned" by the automatic imaging 
process, thus artificially reducing the rms noise of the resulting images. We therefore caution 
that the reader that the noise measured from the images for sources with data only from 
CJF, 141 sources in all, may underestimate the true noise. One may identify these sources 
in Tabled] by their observing dates (i.e., all sources observed before 2004). 

In a few (<1%) cases, the final 5 GHz images were of relatively poor quality and were 
remade within DIFMAP "by hand". The poor quality of the model fits for some of these 
cases was found to be largely due to significant fiux density that was outside of the original 
image in which case larger maps were made. Because the self-calibration process used to 
make the images can produce false point-like sources from peaks in the noise, especially for 
faint objects that are essentially not detected, the maps and visibility data for all sources 
with peak fiux densities <20 mJy beam~^ were inspected, 27 objects in all. The reliability of 
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the maps for these sources was evaluated on a case-by-case basis. The images for a total of 
eleven sources were deemed unreliable and flagged as non-detections and make up about 1% 
of the sample. For all detected sources, the total and peak flux densities at 5 GHz are listed 
in Table [U including those sources observed as part of the CJF and VIPS pilot surveys. 
The 8 sources from the full catalog of 1,127 sources that have or will have data only from 
the MO JAVE survey are not included. The full version of Table [1] is available in electronic 
form only. The version displayed here merely illustrates the general format of the table. 
All Stokes I, Q, and U maps are publicly available in FITS format via the VIPS website, 
http://www.phys.unm.edu/~gbtaylor/VIPS, as well as the VLBA {u,v) data (also in FITS 
format) and PDF and GIF versions of contour plots of the I maps, polarization contours, 
and visibility plots. The full VIPS source catalog (i.e., including the 8 MO JAVE sources 
excluded from Table [T]) is also available via the VIPS website. 



2.3. Source Classification 

The relatively large sample size implies that the VIPS catalog will contain a variety of 
source types. In order to provide some level of information regarding the morphology of the 
VIPS sources, we have divided the sources into four categories (1) point-like sources (PS), 
(2) core-jets, which are subdivided into short jets (SJET) and long jets (LJET), (3) compact 
symmetric object candidates (CSO), and (4) complex sources (CPLX). To deal with the 
relatively large number of sources and to make the classifications as objective as possible, 
we have developed an automated classification procedure based on Gaussian component 
fitting performed within the image plane. The Gaussian fitting was performed within AIPS 
with the SAD task which fits elliptical Gaussians to all sources within an image down to a 
given limiting flux density. Through trial and error, it was found that the best results were 
obtained if the Gaussian fits were performed for sources at ten different flux density limits 
starting with the peak flux density for the image and proceeding down to the 6a level in steps 
evenly spaced in logarithmic flux density. Examples of the Gaussian fits displayed in Fig. 
E] illustrate that the fitting reliably breaks up the objects into distinct components. In Fig. 
m we plot the total flux density contained within the Gaussian components versus the total 
cleaned flux density for all detected sources. While the Gaussian fits clearly overestimate 
the flux density in some instances, in general, the flux density of the Gaussian components 
matches the cleaned flux density relatively well. 

To classify the VIPS sources as objectively as possible, we have developed the following 
algorithm which utilizes the Gaussian components: 

1. If a source has one Gaussian component that contains 95% or more of the total flux 
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density of all its Gaussian components, it is flagged as a single component object. 
Single component objects that are more elongated than the restoring beam used (i.e., 
an axis ratio oih/a < 0.6), are classified as core-jets. Those sources not classified as 
core-jets are classified as point-like (PS). 

2. Sources not flagged as single component sources are flagged as double sources if their 
two brightest Gaussian components contain 95% or more of the total flux density. If the 
flux densities and fitted sizes (assumed to be proportional to ab) of these components 
agree within a factor of 2.5, they are classified as compact symmetric object candidates 
(CSO). If this is not the case, they are classified as core-jets. 

3. Sources not flagged as single or double sources are flagged as multiple component 
sources. The dominant components within each of these sources are identified as the 
brightest sources whose combined flux density is 95% or more of the total flux density. 
For these dominant sources, a line is fit to their relative declinations as a function of 
their relative right ascensions. If the dispersion of the positions of the components 
relative to their center (taken to be the mean position of the components) in the 
direction of this fitted line is a factor of two greater than the dispersion perpendicular 
to the fitted line, the object is classified as a core-jet. Otherwise, the object is classified 
as complex (CPLX). 

4. For all sources classified as core-jets, those longer than 6 mas are classified as long jets 
(LJET). Thos e short er than this limit are classified as short jets (SJET) according to 



Pollack et al. I (120031 ) who found that there may be difference between the polarization 



properties of quasars with jets that are divided into two groups using this limit. 



Following this initial classification, an additional algorithm was run to perform a more de- 
tailed search for CSO candidates with morphologies that are more complex than symmetric 
double sources. This algorithm first identifies sources with two Gaussian components whose 
combined flux density is greater than 80% of the total flux density and whose flux densities 
and fitted sizes (oc ab) agree within a factor of 2.5. This was done to include double sources 
that may have some extended emission that could still be CSOs. Any multiple component 
object classified as LJET with a total length greater than 12 mas whose brightest Gaussian 
component was within 3 mas of the mean position of all components was also reclassified as 
a CSO. This was done to include any source that appears to have a core with significant and 
roughly symmetric diffuse emission on opposite sides of the core oriented along a single axis. 

Following this, the algorithm identifies groups of Gaussian components that overlap on 
the image where in practice, the "groups" are allowed to have as little as one component. 
For this purpose, the outer boundary of each component was defined to be an ellipse with 
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the same position angle as the Gaussian component and major and minor axes equal to 
3a/v81n2 and 36/ v81n2 where a and h and the full widths at half maximum of the Gaussian 
component along the major and minor axes respectively (i.e., effectively 3cr from the center 
of the Gaussian component). Components that had outer boundaries that were defined in 
this manner which intersected were considered to be overlapping components. Using these 
groups, the algorithm identifies the following objects as CSO candidates: (1) objects with 
two groups that contain 80% of the total flux density and whose flux densities agree within a 
factor of 2.5, (2) objects that have at least two multiple component groups where the group 
closest to the center (equal to the mean component position) has only one component (i.e., it 
is likely the core at the center of the CSO), and (3) objects with more than two groups where 
the brightest component is closest to the center. Images of examples of sources reclassified 
as CSO candidates by this algorithm are displayed in Fig. [71 

To test the quality of the automatic classifications, visual inspection of the I image for 
each source was performed. The "by-eye" classifications derived from these visual inspections 
agreed with the automatic classifications in 99% of the cases for both point-like objects 
and short jets, in 93% of the cases for long jets, in 87% of the cases for CSO candidates, 
and in 71% of the cases for complex objects. The larger discrepancy between the two 
classifications for CSO candidates and complex objects appears to have more to do with 
peak flux density than with morphology. In general, the by-eye classifications agree with the 
automatic classifications more frequently for brighter objects. This is illustrated in Fig. [5] 
where we have plotted the fraction of sources for which the two classifications agree within 
bins of peak 5 GHz flux density. These results show that for sources with peak flux densities 
greater than ~60 mJy, the two classifications agree in 95% of the cases. The median peak 
flux densities for the CPLX, CSO, LJET, SJET, and PS sources are 24, 54, 84, 92, and 
102 mJy beam~^ respectively, implying that higher fractions of CSO and CPLX sources are 
misclassified due to their relatively low flux densities and not their complex morphologies. 

Based on the comparison with the by-eye classifications, the automatic classification 
algorithm appears to provide reliable and objective source types and may be used successfully 
with follow-up VLBI observations of the VIPS sample or with other VLBI imaging surveys. 
However, since the performance of the algorithm is lower for the typically fainter CSO 
candidates and complex sources which are among the rarest and most interesting sources in 
the sample, we have elected to make our by-eye classifications available as well to facilitate 
follow-up observations of these types of sources. To this end, we have listed the automatically 
determined source types in Table [T] along with the by-eye classifications for those sources 
where the two classifications disagreed. For objects flagged as non-detections, the type is 
listed as ND. A summary of the number of sources in each category as well as the fraction of 
sources of each type with detected polarized flux density is contained in Table [2l For Table 
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ls] and the remainder of the paper, the source types used include the by-eye re-classifications 
where applicable. 

For each object, the number of dominant Gaussian components is listed in Table [H For 
each object, these dominant Gaussian components were also used to measure a radius equal 
to -\/a6 for single component sources and equal to the average distance of the components 
from their centroid for double and multiple component sources which is listed in Table [H For 
double and multiple sources, the maximum separation among the components of each object 
is also listed in Table [H along with a position angle for any objects classified as core-jets 
computed using a linear fit to the relative right ascensions and declinations of the dominant 
Gaussian components. 

2.4. Polarization Properties 

To exploit the relatively large number of source with full polarization data provided by 
VIPS, we have developed the following automated method for measuring the polarization 
properties of different components of each source. First, for each source with detected 
(> 5a) polarized flux density, we constructed three image masks, one using all the Gaussian 
components flt to the I image (see §2.3), one using only the brightest Gaussian component, 
and one that is the difference between these two masks. We then constructed another image 
mask using the polarized intensity and noise images produced by the AIPS task COMB 
using the Q and U images from the 4994 plus 5095 MHz pair (see §2.2) including error 
biasing. This mask was made by setting pixels with signal-to-noise ratios >5 to unity and 
the remaining pixels to zero. We then constructed three composite masks by multiplying 
the polarized intensity mask by the three masks made using the Gaussian components and 
used them with the polarized flux density and I images to measure the mean polarization 
fraction, fpoi, or the ratio of polarized to total intensity. Using the three separate masks, 
we obtained measurements of fpoi for the entire object, the object's core, which we assumed 
to be represented by the brightest Gaussian component, and the regions outside the core. 
In the majority (~55%) of sources with more than one dominant Gaussian component, 
the polarized flux density is only found within the core. About 35% of these sources have 
detected polarized flux density both within and outside the core while about ~10% have 
detected polarized intensity only in the regions outside the core. 

We also used the same three composite masks to measure a polarized intensity-weighted 
mean EVPA for the whole object, the core, and the regions outside the core using the po- 
larized flux density image and the polarization angle image, also produced using the COMB 
procedure. For sources with both positive and negative EVPA values in their polarization 
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angle images, care was taken to ensure that the mean EVPA was computed properly to 
ensure that, for instance, for a source with EVPA values near both 90° and -90°, the mean 
EVPA was near either 90° or -90° and not 0°. Specifically, in each of these instances, mean 
EVPA values were computed separately for the positive and negative pixels on the polar- 
ization angle image. If the difference between these two mean values was less than 90°, a 
polarized intensity-weighted EVPA was computed using all the pixels from the polarization 
angle image. If the two mean values differed by more than 90°, 180° was added to each 
negative pixel on the polarized image after which the polarized intensity-weight mean EVPA 
was computed. If this mean EVPA was greater than 90°, 180° was subtracted from it so that 
all of the polarized intensity- weighted EVPA values would be between —90° and 90°. All 
polarization related quantities are listed in Table [3] for all objects with detected polarized 
flux density. As with Table [U we have included only the first twenty sources here to provide 
an example of the table format, and the entire table is available in electronic form only. We 
have used the flux densities from the I, Q, U, and polarized intensity images along with the 
rms values measured from the I, Q, and U images to estimate the uncertainties in both fpoi 
and X ci-nd find that the typical errors in these quantities are 0.003 and 3° respectively. 



3. Results and Future Work 

The compilation of images and polarization data that make up VIPS constitute the 
largest such database of AGN to date. A number of scientific endeavors to explore the 
nature of AGN are possible with these data. Here, we will briefly explore the most basic 
properties of the sample and will leave more detailed analysis for subsequent papers. 



3.1. Comparison with Optical Data 



Since the VIPS sample was chosen to lie on the SDSS survey footprint, there are existing 
optical data for nearly all of our sources. Throu gh DR5, 997 (88%) hav e optical magnitudes 
measured by the SDSS photometric pipeline (see IStoughton et al. II2002I ) and 356 (32%) have 
SDSS spectra. In Fig. [SI we have plotted the 5 GHz flux densities measured from the VIPS 
images versus the SDSS i-band magnitudes separately for each of the five source categories. 
For those VIPS sources with SPSS sp ectra, we have used the sources' redshifts and the 
K-corrections of [Richards et al. I (120061 ) to correct their z-band magnitudes to a redshift of 
z=0 and have plotted these sources separately in the right panels of Fig. [SI For both the 
observed and K-corrected i-band magnitudes, there is no clear trend between the 5 GHz 
flux densities and optical magnitudes for any of the five source categories. This illustrates 
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the need for further optical follow-up observations to obtain redshifts for as many VIPS 
sources as possible so that consistent distance measurements may be obtained without any 
bias introduced by the selection function of the SDSS and any additional optical selection 
effects. Optical spectra are currently being obtained for all VIPS sources not targeted for 
spectroscopy by the SDSS. A detailed discussion of the results of this work will be presented 
in a subsequent paper. 

One of the most useful properties of broad-line AGN that can be estimated using optical 
spectra is the virial mass of the central black hole, assumed to be equal to G~^RblbV^ where 
Rblr is the radius of the broad line region and V is th e velocity width at half maximum 



of the broad optical emission lines. iKaspi et al. I (120001 ) have demonstrated that Rblr is 



strongly correlated with the monochromatic continuum luminosity at 5100 A, L^iqq. Using 
their observed correlation, one can use the velocity width of the H/3 emission line, V{H(3), 
with L5100 to estimate the central black hole mass, Mbh- With such a large sample of high 
resolution radio frequency images of sources that have or will have optical spectra, we are 
in a good position to explore any relation between Mbh and the radio frequency luminosity 
of the cores of AGN or other properties. 

Using th e Gaussian fits to the emission lines performed by the SDSS spectroscopic 



pipeline (see IStoughton et al. II2002I ). we have computed V{H/3) for all sources that have 
> 3cr detections of the B.f5 emission line and which are at low enough redshifts that values 
for rest-frame L5100 could be determined, 62 sources in all. In Fig. [HI we have plotted the 
estimated values of Mbh for these 62 sources versus their total 5 GHz luminosity densities and 
versus the luminosity density of the core component of each point-like and core-jet source 
(57 sources in all), which we take to be the luminosity density of the brightest Gaussian 
component. In both cases, we have included errors in the luminosity densities that reflect 
the range in rest frame 5 GHz luminosity density expected for power-law spectra with slopes 
ranging from —0.5 to 0.5. We find a slight correlation between black hole mass and 5 GHz 
luminosity density with a Spearman rank-order correlation coefficient of 0.4. The probability 
of getting this result by chance, however, is about 50%, indicating that the trend is weak at 
best. The results are nearly the same if we only consider sources with V{H(3) >2,000 km 
s~^, or if we only consider the core luminosity for point-like and core-jet sources. However, 
with the completion of the optical follow-up, we will be able to explore this issue much more 
thoroughly using velocity widths and radio luminosities for ~3 times as many sources. 
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3.2. Polarization Results 

3.2.1. Fractional Polarization 

Perhaps the most distinguishing aspect of VIPS is the number of sources with detected 
polarized flux density, 393 sources in all and about 37% of all newly observed sources (i.e., 
excluding sources with data only from the CJF or VIPS pilot surveys). This puts us in 
a position to be able to accurately measure, among a relatively large sample of AGN, the 
distribution of fractional polarization, /poj, for different source components (see §2.4) and 
different source types. In Fig. [TOl we have plotted the fpoi distributions for all newly observed 
sources for the entire objects, the objects' cores, and the regions outside the cores. Using 
all regions of polarized intensity from each object, the median value of fpoi is about 5% 
and ranges from 1-20%. The results are similar if only the cores are used. However, f^oi is 
significantly larger on average for regions outside the cores with a median value of ~17% and 
a range of about 7-100%. These regions are generally the jets of core-jet systems since 90% 
of systems with polarized intensity detected outside the cores are classified as core-jets with 
the remaining 10% being CSO candidates and complex sources. This implies that relatively 
speaking, jets are more strongly polarized than cores on average within core-jet systems. 
In fact, out of the 92 sources with polarized flux density detected both within and outside 
the cores, only one source has a value of fpoi that is larger for the core than for the regions 
outside the core. 

To explore any trend between fpoi and source type, we have displayed a so called "box- 
and-whisker" plot of fpoi versus source type in the lower panel of Fig. [TD], excluding CPLX 
sources for which there were only two sources with detected polarized flux density. These 
plots indicate that the median, upper and lower quartiles, and extreme values of fpoi are all 
larger for LJET sources than for both point-like and short jet sources, regardless of whether 
the entire object, the core, or the regions outside the core are used to compute fpoi- The 
median values for CSO sources are similar to those for LJET sources. One should keep in 
mind, however, that objects classified as CSO are CSO candidates and that many (if not all) 
of these candidates with detected polarized flux density may in fact be core-jet systems. The 
results in Table [2] imply that in fact, detected polarized intensity is relatively rare among 
the CSO candidate sources with only 15 out of 103 sources having significant polarized flux 
density. By comparison, about 40% of LJET sources have detected polarized flux density, 
the highest of any source type, which is at least in part due to the relatively large values of 
fpoi found for these sources. 
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3.2.2. EVP A and Jet Direction 

To explore the possibility of a relationship between the direction of core magnetic fields 
and the direction of the cores' associated jets, we have computed the absolute difference 
between the polarization angle and jet position angle of each source with detected polarized 
flux density classified as a core-jet (see Fig. [TT] for some examples of polarized sources with 
jets). For each of these sources, we take the measurement of the EVPA, x, outside the 
core (see §2.4) to be the EVPA of the jet component, or Xjet- For this analysis, we have 
excluded any source with a single dominant Gaussian component that was classified as a 
core-jet because it was more elongated than the restoring beam used. This was done for two 
reasons. First, while such objects are too elongated to be unresolved point sources, they 
have not been resolved into separate components and measuring the polarization properties 
of their core and jet components separately would be difficult. Second, since these sources 
are only marginally resolved, their jet position angles will tend to the position angle of the 
restoring beam (i.e., PAjet = 0). This only effects SJET sources, slightly less than half of 
which each have a single dominant Gaussian component. 

In Fig. [12], we have plotted the distributions for |x — PAjet|, \Xcore — PAjet|, and \xjet — 
PAjetl for all core-jet systems and separately for SJET and LJET sources. For each of these 
distributions, we have used a K-S test to compute the probability that the distribution was 
drawn from a flat distribution and have printed the results in the corresponding panels of Fig. 
[T2]as Pfiat- For all core-jet and LJET sources, the |x — PAjetl and \xcore — 'P^jet\ distributions 
are noticeably peaked near 90°; the probability that each was drawn from a fiat distribution 
is <0.001%. This implies that these peaks are likely the result of a real tendency for the core 
EVPAs and jet position angles within core-jet systems t o be perpendicular t o one another 



which is consistent with what was found for quasars by iPoUack et al. I (120031 ). The results 
for the remaining distributions are more marginal but hint that the same is true for SJET 
sources and that jet position angles may also tend to be anti-aligned with jet EVPAs, but 
much less frequently than with core EVPAs. 

To explore the tendency for Xcore and PAjg* to be anti-ahgned implied by the distribution 
for all core-jet and LJET sources shown in Fig. [121 we have constructed the following simple 
model. First, we assume that a significant fraction of polarized core-jet systems have intrinsic 
polarized flux densities, Pq, with components Qq and Uq which have been altered by a 
combination of observational errors and Faraday rotation and which have EVPAs that are 
perpendicular to the jet axes. We then define a new coordinate system within the (Q,U) 
plane by rotating the Q and U coordinates such that the new coordinates, Q' and U' , are 
perpendicular and parallel to the jet axis respectively, i.e., Qo = and Uq = Pq. To simulate 
the influences of observational error and Faraday rotation, we assumed that magnitude of 
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both these effects is the same for the Q' and V components of the polarized flux density. We 
then assumed that these two effects are additive and that the distribution of the appropriate 
additive factors can be approximated by a single Gaussian function such that the Q' and V 
components are given by 

Q' = Rgqct (2) 
f/' = Po + Rguct (3) 

where Rgq and Rqu are two separate random numbers drawn from unit normal distribu- 
tions and a represents the rms uncertainty in both the Q' and U' flux densities caused by 
observational errors and Faraday rotation. For convenience, we also define a parameter 
c = P^/a so that we may use two sets of random numbers drawn from unit normal distribu- 
tions to compute model distributions for \xcore — PAjet|± = tan~^(t/'/(5')/2. We also allow 
for the possibility that some fraction of core jet systems have core EVPAs and jet position 
angles that are completely unrelated by computing a second model distribution according 
to \Xcore — ^^jet\random = ^^^yT^ {Rqu I Rgq) / ■ Usiug 10^ pairs of raudom numbers, we used 
these two model distributions to iteratively solve for the best fitting values of c and the 
fraction of sources with anti-aligned EVPAs and jet position angles. The observed and best 
fitting distributions are plotted in Fig. [T31 the best fit model is for c = 3.0 and implies that 
the fraction of core-jets with anti-aligned EVPAs and jet position angles is 0.24. An estimate 
of the covariance matrix for these parameters was computed, and it was found that the two 
parameters are moderately anti-correlated with a correlation coefficient of —0.66. 

The best fitting value of c implies that on average, ct/Pq —0.33. For the newly observed 
VIPS core-jet systems with detected polarized flux density, the median ratio of the rms errors 
measured from the Q and U images to the peak polarized flux density is ~0.1. This implies 
that in order for the model results to be reasonable, Faraday rotation must dominate the 
uncertainty in the EVPAs, contributing nearly 80% of the uncertainty in the Q sand U flux 
densities. This is not only reasonable, but expe cted since the RMs of quasars have been 



observed to be ~500 to a few thousand rad m~^ fjZavala &: Taylor II2004I ). corresponding to 
rotations of at least ~100° at a frequency of 5 GHz. A proper and more detailed exploration of 
the relationship between core EVPA and jet position angle will require follow-up observations 
at multiple frequencies so that RM values may be obtained for our core-jet sources and the 
appropriate corrections for Faraday rotation can be made. 



3.3. Additional Follow-up 



With this initial data set, we now have the ability to design follow-up experiments to 
yield even more information about the nature of AGN. In particular, with a relatively large 
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sample of good CSO candidates, we are in a good position to use follow-up VLBI observa- 
tions at other frequencies to confirm whether or not these objects are in fact CSOs and to 
use multiple epochs to explore how these objects evolve. Among the CSO candidates and 
complex sources, we have identified ~20 sources that are good candidates for small separa- 



tion s upermassive binary black hole (SBBH) systems similar to 0402-1-379 (IRodriguez et al. 



20061 ). A good example of one of these sources, J10019-I-5540, can be seen in the middle 



panels of Fig. [31 Ongoing VLBA follow-up observations at 5, 8, and 15 GHz will allow us 
to confirm whether these systems are indeed compact SBBH systems. We will also be able 
to obtain rotation measures for all polarized core-jets (~300 sources) using follow-up VLBI 
observations at additional frequencies, allowing for a much more thorough exploration of the 
relation between magnetic fields and jet activity. Currently, follow-up observations of all 
core-jets with detected polarized flux density, all CSO candidates, and all complex sources 
not included in the ongoing SBBH candidate follow-up program are planned with the VLBA 
at 5, 8, and 15 GHz. Finally, future observations of core-jet sources with GLAST will allow 
us to constrain the physics involved with jets even further. 
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Table 1. Source Properties 
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Note. — CoL (1): VIPS source number. CoL (2): VIPS source name. CoL (3): Right ascension (J2000). Col. (4): Declination (J2000). Col. (5): UT date of the 
observations. Col. (6): The flux density at 8.5 GHz from the CLASS survey. Col. (7): The total cleaned flux density from the 5 GHz VLBA map. Col. (8): The 
pealc flux density from the 5 GHz VLBA map. Col. (9): The rms noise of the 5 GHz VLBA image. Col. (10): The number of dominant Gaussian components (i.e., 
that contain more than 95% of the total flux) fit to the 5 GHz VLBA map (see §2.3). Col. (11): The source type derived from the automated Gaussian component 

classification (sec §2.3). Col. (12): The source type determined by visual inspection of the I image for sources where the "by-eye" and automatic classifications 
disagree. Col. (13): The mean radius (i.e., mean distances from the mean component position) for the ensamble of dominant Gaussian components. Col. (14): The 
maximum separation among the dominant Gaussian components. Col. (15): The jet position angle (measured from north through east). 
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Table 2. Summary of Source Types 
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Table 3. Source Properties 
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64 
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0.033 
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Note. — Col. (1): VIPS source number. Col. (2): VIPS source name. Col. 
(3): Average 5 GHz fractional polarization (see §2.4). Col. (4): Average frac- 
tional polarization within the brightest Gaussian component. Col. (5): Average 
fractional polarization outside the brightest Gaussian component. Col. (6): Po- 
larized intensity-weighted mean electric vector position angle (EVPA) at 5 GHz 
(see §2.4). Col. (7): Polarized intensity-weighted mean EVPA within the bright- 
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est Gaussian component. Col. (8): Polarized intensity- weighted mean EVPA 
outside the brightest Gaussian component. 




Fig. 1. — The sky coverage of the VIPS source catalog superimposed on the sky coverage 
of the original SDSS footprint. The light gray shaded area represents the original SDSS 
footprint; the dark gray shaded are represents the area with SDSS images through the 
fifth data release (DR5); the white open circles represent the spectroscopic plates observed 
through DR5. The solid black points represent VIPS sources; the open black points represent 
those VIPS sources that are outside the SDSS DR5 imaging area. 



- 25 - 




Fig. 2. — The ratio of the rms noise measured from the 5 GHz image, rmsimage; to the 
theoretical noise computed according to equation (1), rmstheory, for VIPS sources observed 
in 2006 (sohd points ) and VIPS sources imaged using data from the CJF survey (open points, 



Britzen et al. 



20031 ) versus the peak signal-to-noise ratio from the 5 GHz image. The median 



ratios of rmsimage/rmstheory are represented by a dashed line for the VIPS sources observed 
in 2006 and by a dotted line for CJF/VIPS sources. 
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Fig. 3. — For three VIPS sources, the 5 GHz I maps (left column), the Gaussian fits to the I 
maps (middle column), and the residuals for the Gaussian component fits (right column) for 
a core-jet system (upper panels), a complex system (middle panels), and a CSO candidate 
(lower panels). 
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Fig. 4. — The total flux density at 5 GHz contained within the components from the Gaussian 
fits (see §2.2) versus the total cleaned flux density. The dashed line is not a fit, but simply 
the line expected for Fqauss—Fclean plotted for reference. 
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Fig. 5. — Within bins of peak 5 GHz flux density, the fraction of sources for which the 
"bye-eye" classification agrees with the automatic classification algorithm (see §2.3). The 
median peak 5 GHz flux densities for the flve sources categories are displayed as vertical 
dotted lines. 
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Fig. 6. — An example of a point-like object (upper left), a single component core-jet (upper 
right), a double component core-jet (middle left), a multiple component long jet (middle 
right), a CSO candidate (bottom left), and a complex source (bottom right) as determined 
by the automated classification algorithm detailed in §2.2. Here, "sigma" and "peak" refer 
to the reduced between the CLEAN components and the data and the peak cleaned flux 
density in Jy beam"-^ respectively. 
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Fig. 7. — Examples of four sources reclassified as CSO candidates by the additional special- 
ized CSO classification algorithm (see §2.3). Here, "sigma" and "peak" refer to the reduced 
between the CLEAN components and the data and the peak cleaned flux density in Jy 
beam~^ respectively. 
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Fig. 8. — The 5 GHz flux density measured from the VIPS images versus the SDSS i-band 
magnitude for each of the flve source categories (see §2.3). In the left panels, the observed 
2-band SDSS magnitudes are used; in the right panels, only the sources with SDSS spectra 
ar e included and their i-b and magnitudes have been K-corrected to z=0 using the corrections 
of [Richards et al. I J2006h . 
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Fig. 9. — For sources with SDSS spectra and > 3a detections of the H/? emission line, the 
estimated viral mass of the central black hole (sec §3) versus the total 5 GHz luminosity 
density (upper) and the luminosity density of the core component for point-like and core-jet 
sources (lower) taken to be the luminosity density of the brightest Gaussian component (see 
§2.3). In both panels, the error in each luminosity density reflects the range in rest-frame 5 
GHz luminosity densities expected for a power-law slope of —0.5 < a < 0.5. The Spearman 
rank-order correlation coefficient between Mbh and L(5 GHz) is 0.4. The probability of 
getting this result by chance is ~50%, implying that the observed trend is rather weak. 
Similar results were obtained when only sources with H/3 velocity widths, V{Hf3), greater 
than 2,000 km s~^ (represented as stars in both panels) were considered and when only 
Lcore{5 GHz) was considered for point-like and core-jet sources. 
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Fig. 10. — Upper: The distribution of the mean 5 GHz polarization fraction, fpoi, measured 
for the entire object (solid line), the core (i.e., the brightest Gaussian component; dashed 
line), and the regions outside the core (dotted line; see §2.4). Lower: a so called "box-and- 
whisker" representation of the three fpoi distributions for each source category (see §2.3) 
excluding complex sources for which there were only two sources with detected polarized 
intensity. The lower and upper boundaries of each box represent the 25*^ and 75*^ percentiles 
respectively. The horizontal line within each box represents the median, and the dashed lines 
extend to the extrema of each distribution. 
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Fig. 11. — Examples of six objects identified as core-jets (see §2.3) that have detected 
polarized fiux density. In each plot, the contours correspond to the total intensity and the 
orientation of the lines is taken from the polarization angle image (see §2.4) and represent 
the EVPA without any correction for Faraday rotation. The length of each line represents 
the ratio of the polarized to total intensity; the scale for these lines is listed below each plot. 
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Fig. 12. — The distributions for the absolute difference between the jet position angle and the 
polarized intensity-weighted mean EVPA, x, measured for each object (left panels), the core 
of each object (middle panels), and the jet component of each object which we take to be the 
regions outside the core (right panels; see §2.4). Distributions are displayed separately for 
all core-jet systems (upper panels), SJET sources (middle panels), and LJET sources (lower 
panels) . The probability that each distribution was drawn from a flat distribution computed 
using a K-S test is displayed in the appropriate panel as Pfiat', values of Pfiat < 0.001% are 
listed as Pfiat — 0. 
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Fig. 13. — The observed (solid line and error bars) and best fitting model (dotted line; see 
§3.2) distributions for the absolute difference between the core EVPA and the jet position 
angle for core-jet systems with detected polarized flux density. The best fitting model implies 
that about 24% of core-jet systems have EVPAs that are perpendicular to their jet position 
angles and that the ratio of the uncertainties in the Q and U flux densities for these systems 
are dominated by Faraday rotation (see §3.2). 



